Doped coupled frustrated spin-| chains with four-spin exchange 
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The role of various magnetic inter-chain couplings is investigated by numerical methods in doped 
frustrated quantum spin chains. A non-magnetic dopant introduced in a gapped spin chain releases 
a free spin-i soliton. The formation of a local magnetic moment is analyzed in term of soliton con- 
finement. A four-spin coupling which might originate from cyclic exchange is shown to produce such 
a confinement. Dopants on different chains experience an effective space-extended non-frustrating 
pairwise spin interaction. 
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Strong interest for quasi-one-dimensional frustrated 
spin chains has been revived since the discovery in 
1993 of the first inorganic spin-Peierls (SP) compound 
CuGeC>3 which exhibits a transition Q] towards a gapped 
dimerised ground state (GS). Some vanadate oxides such 
as LiV2 05 are also excellent realizations of weakly inter- 
acting frustrated spin- 1/2 chains Whether magneto- 
elastic couplings play a role at low temperature in 
L1V2O5 is not vet clear although it has been investigated 
theoretically pj- Replacing a spin-^ in a spontaneously 
dimerised (isolated) spin chain by a non magnetic dopant 
(described as an inert site) liberates a free spin i , named 



a soliton, which do not bind to the dopant 4]. The soli- 
ton can be depicted as a single unpaired spin (domain) 
separating two dimer configurations 0] . The physical pic- 
ture is completely different when a static bond dimerisa- 
tion exists and produces an attractive potential between 
the soliton and the dopant |4],|5( and consequently leads, 
under doping, to the formation of local magnetic mo- 
ments 0, as well as a rapid suppression of the spin 
gap . However, a coupling to a purely one-dimensional 
(ID) adiabatic lattice |8( does not produce confinement 
in contrast to more realistic models including an elastic 
inter-chain coupling (to mimic 2D or 3D lattices) 0, 0] . 
Doping quasi-lD dimerised spin chains is realized ex- 
perimentally by substituting a small fraction of Copper 
atoms by Zinc (or Magnesium) atoms in CuGeC>3 jlfj. 



The motivation of this Letter is to investigate by nu- 
merical methods based on Exact Diagonalisations (ED) 
of finite chains (see Ref. |9j) the role of various inter-chain 
magnetic couplings on soliton confinement. Such cou- 
plings are expected to become crucial in the case where 
magneto-elastic effects remain small. Indeed, in that 
case, we argue that the formation of magnetic moments 
in doped gapped frustrated spin chains is then primarily 
controlled by 4-spin interactions. 



Let us first consider a model of coupled frustrated spin- 
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FIG. 1: Schematic picture of the coupled chains model with 
nearest neighbor, next-nearest neighbor, inter-chain and 4- 
spins couplings Ji, J2 = aJi, Jj_, and J4. Full (resp. open) 
circles stand for spin-i sites (resp. non-magnetic dopants). 



h antiferromagnetic (AF) chains, (see FigGJ 

i.a 

], (i) 

where i is a lattice index along the chain of size L and a 
labels the M chains (L and M chosen to be even). Peri- 
odic boundary conditions will be assumed in both direc- 
tions. In the absence of inter-chain coupling ( Jj_ = 0) the 
behavior of the un-doped chain is controlled by the frus- 
tration a. From a — to a c — 0.241167 [ill], the chain 
is in a critical regime with a gapless excitation spectrum. 
Larger frustration stabilizes a spontaneously dimerised 
phase characterized by a gap opening. At a = 0.5, the 
so-called Majumdar-Ghosh |T3 (MG) point, the 2-fold 
degenerate Ground State (GS) consists in the product 
of singlets located either on odd or on even bonds. Be- 
yond the MG point, the short-range correlations become 
incommensurate . 

The inter-chain AF coupling J± can stabilize a 
Neel ordered phase under some conditions. Following 
Schulz 0|, a mean-field (MF) treatment of the inter- 
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FIG. 2: SP-AF phase diagram in the (a, J±) plane from ED 
of chains of length up to 16 sites. Symbols correspond to dif- 
ferent values of J 4 > as indicated on plot. Typically, FSE 
are smaller than the size of the symbols. The computed tran- 
sition lines are extended by tentative transition lines (dashed 
lines) in the region where FSE become large. At J4 = we 
have plotted a few points in the vicinity of the MG point for 
L — 12 and L — 16. 
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FIG. 3: Local magnetization Mj E for L x M=16 x 8 coupled 
chains with one dopant D (shown by arrow) located at a = 1, 
i = 16 in the dimerised phase (a = 0.5, J± = 0.1). Circles 
correspond to J4 = (shown up to the third neighboring 
chain of the doped one) and squares (crosses) to J4 = 0.01 
(J4 = 0.08). The coupling across the dopant has been set to 
for convenience. 



chain coupling is performed, 

H e s((X, Jj_) = ^^[5i,a • Si + i. a + aSi ia ■ Oj+2, a 
i,a 

+ hi,aSi a — J±(Si a )(S* a+1 )} 7 (2) 

with 

h i , a = Jx((Sl a+1 ) + (Sl a _ 1 )), (3) 

the local magnetic field computed self-consistently. In 
the absence of dopant, we expect an homogeneous AF 
phase characterized by a self-consistent staggered mag- 
netization (Sf a ) — (—iy +a m so that the system re- 
duces to a single chain under a staggered magnetic field 
hi = ±2(—iyj±m. In the absence of frustration (a = 0), 
it was shown that m ~ \fJZ By solving the self- 

consistency condition using ED of finite chains with up to 
16 sites (supplemented by a finite size scaling analysis) 
the transition line J± = (a) (see Fig. |2J separating 
the dimerised SP phase (m = 0) and the AF ordered 
phase (for which m ^ 0) has been obtained in agreement 
with field theoretic approaches 0|. Finite size effects 
are small in the gapped regime and especially at the MG 
point (Jj_(ot = 0.5) ~ 0.11 for all sizes). Note also that 
numerical data suggest that the AF order sets up at ar- 
bitrary small coupling when a < a c with a clear finite 
size scaling J]_(L) cx 1/L at small a. 

A non-magnetic dopant is described here as an inert 
site decoupled from its neighbors. Under doping the sys- 
tem becomes non-homogeneous so that we define a local 
mean staggered magnetization, 

AC S = \(-lT +a (2(Sl a ) - (S? +1<a ) {SU, a ))- (4) 



Following the method used in Ref. |2( , the MF equations 
are solved self-consistently on finite L x M clusters and 
lead to a non-uniform solution. At each step of the MF it- 
eration procedure, we use Lanczos ED techniques to treat 
exactly (although independently) the M non- equivalent 
finite chains and compute (S* a ) for the next iteration 
step until convergence is eventually achieved. 

We first consider the case of a single dopant. Whereas 
in the AF phase the magnetization profile is weakly en- 
hanced in the vicinity of the dopant 0, Ht| , in the 
SP gapped phase the soliton remains de-confined as can 
be seen from FigG3 Note that the inter-chain coupling 
induces a " polarization cloud" with strong antiferromag- 
netic correlations in the neighbor chains of the doped 
one. 

The fact that our treatment of J± does not produce 
confinement is clear at the MG point where the dopant 
site separates two regions with out-of-phase dimer cover- 
age. In such a state, the local magnetization (Sf) iden- 
tically vanishes. Hence, we expect J±Sf a (S£ a ±i) to be 
irrelevant. On the contrary, terms like (Si ■ Si+i) are or 
—3/4 in the MG state and, hence, are relevant. There- 
fore, a 4-spin coupling Ji(S ita ■ S i+ i, a )(Si±i iQ+ i • Si, a +i) 
which originates from cyclic exchange |18j is of crucial 
importance As done for J± we have applied a 

mean-field treatment to the J4 coupling leading to a (self- 
consistent) modulation of the intra-chain nearest neigh- 
bor couplings of the Halmiltonian, 

H e g(a, J±, J4) = y^[(l + SJi ya )Si : a ■ Si+l.a 
i,a 

+aSi, a ■ §i+2, a + hi ia S? a + constant , (5) 
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with hi <a given by Eq.@ and 

8Ji,a = J4,{{Si,a+l ' Si+l.a+l) + (S^a-l ' Sj+l )0 -l)}- (6) 

As above, the numerical treatment again consists in a 
recursive diagonalisation of the M spin chains in the 
presence of the (non-uniform) MF magnetic fields and 
exchange modulations produced by its neighbors. Note 
that, as before, the effective coupling between the chains 
appears only through the MF self-consistency conditions 
so that, in practice, each chain Hamiltonian can be Lanc- 
zos diagonalised independently at each iteration. 

Let us first consider the relative stability of the AF 
phase with respect to the dimerised phase in the un- 
doped system with 4-spin coupling. Since all chains be- 
come equivalent we have to deal with an effective single 
dimerised frustrated chain problem in a staggered mag- 
netic field. If the chains are not dimerised, 5Ji^ a is con- 
stant and the exchange is just renormalised. On the other 
hand, if the SP phase is stable, each chain displays the 
same dimerised pattern when J4 < whereas dimer or- 
der is staggered in the transverse direction when J4 > 0. 
Note that, apart from special features (see below), physi- 
cal properties for positive or negative J a are quite similar. 
Note also that the dimerised GS would be 2 M -fold degen- 
erate when J4 = (each chain is independently two-fold 
degenerate) while the degeneracy is reduced to 2 when 
J4 7^ 0. When a = J± = Q but J 4 7^ each chain spon- 
taneously dimerises and a gap opens up. Consequently, 
a finite value of the AF inter-chain coupling is necessary 
to drive the system into the AF ordered phase [2o|. The 
frustration a stabilizes further the dimerised phase with 
respect to the AF one, the critical J±(a) increasing with 
increasing a as seen from the phase diagram shown in 
Fig. |21 Note that finite size effects (FSE) become large 
at small Ja and a and accurate extrapolations are not 
possible there. 

Contrary to the previous case the spin-^ soliton re- 
leased in a dimerised state by doping experiences a con- 
fining string which binds it to the dopant. Indeed, the 
J4 interaction stabilizes (if one assumes an infinite num- 
ber of chains) a given dimerisation pattern (over the two 
possible). Fig.[3]shows the enhancement of the staggered 
magnetization J3J in the doped chain in the vicinity of the 
dopant when J4 7^ 0. Note that the dopant side where 
the soliton is bound is imposed by the (arbitrary) sign of 
the bulk dimerisation. A local spin-1/2 magnetization is 
then expected in the vicinity of the dopant. 

In order to "measure" the strength of confinement, 
we define an averaged soliton-dopant distance as £ = 
Ei^K^I/Ei \{S*)\ so that £ = L/2 in the absence of 
confinement (J4 = 0) and £ converges to a finite value 
otherwise when L — ► 00. On FigQ] £ is plotted ver- 
sus J4 for 2 different system lengths and a — 0.5 and 
J_l = 0.1. FSE decrease for increasing J4. Note that 
£(^4) 7^ £( — J a) an d a power law []| with different expo- 
nents rj is expected when J4 — > 0. A fit gives 77 ~ 0.33 
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FIG. 4: ED data of the soliton average position vs J4 calcu- 
lated for a = 0.5 and J± — 0.1. Different symbols are used 
for L x M = 12 x 6 and 16 x 8 clusters. The long-dashed line 
is a power-law fit (see text). Inset shows the magnetization 
profile in the doped (a = 1) chain at J4 = 0.08, ie £ ~ 2.5. 



if Ja < and r\ — 0.50 for J 4 > (FigQJ. This asym- 
metry can be understood from opposite renormalisations 
of Ji for different signs of J a- Indeed, if J a < then 
8Ji : a > and the nearest neighbor MF exchange becomes 
larger than the bare one. Opposite effects are induced by 
Ja > 0. 

We now turn to the investigation of the effective in- 
teraction between dopants. We consider here a system 
of coupled chains with two dopants (see Fig. Q}. Each 
dopant releases an effective spin | which is localized at 
a distance ~ £ from it due to the confining potential set 
by Ja- We define an effective pairwise interaction J eff as 
the energy difference of the S — 1 and the S = GS. 
When J cff = E(S = 1) - E(S = 0) is positive (negative) 
the spin interaction is AF (ferromagnetic). Let us first 
consider the case of two dopants in the same chain, (i) 
When the two vacancies are on the same sub-lattice the 
moments experience a very small ferromagnetic J off < 
as seen in Fig. [5f A) so that the two effective spins \ 
are almost free, (ii) When the two vacancies sit on dif- 
ferent sub-lattices, Ai is odd and the effective coupling 
is AF with a magnitude close to the singlet-triplet gap. 
Fig- EtA) shows that the decay of J eff with distance is in 
fact very slow for such a configuration. Physically, this 
result shows that a soliton and an anti-soliton on the 
same chain and different sublattices tend to recombine. 

The behavior of the pairwise interaction of two dopants 
located on different chains (Aa = 1,2) is shown on 
Fig- 13 B) & (C) for J 4 > 0. When dopants are on oppo- 
site sub-lattices the effective interaction is antiferromag- 
netic. At small dopant separation J cS (Ai) increases with 
the dopant separation as the overlap between the two AF 
clouds increases until Ai ~ 2£. For larger separation, 
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FIG. 5: Magnitude of the effective magnetic coupling between 
two impurities located either on the same chain (A) (Aa — 0) 
or on different ones (B) & (C) (Aa = 1,2) vs the dopant 
separation Ai in a system of size L x M = 16 x 8 with a — 
0.5, J± = 0.1, and J4 = 0.08. Closed (resp. open) symbols 
correspond to AF (F) interactions. 
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FIG. 6: Local staggered magnetization ^ g for 2 dopants 
(-Di & -D2) on neighboring chains. The labels (a)-(d) refer to 
the data points shown on Fig. EI-B) corresponding to different 
relative positions of the dopants. Arrows stand for the dopant 
positions in the chain direction and full (open) symbols cor- 
respond to the chain doped with D2 (Di). 



was found in second order perturbation in J± \2c 



J cS (Ai) decays rapidly. Note that the released spin-i 
solitons bind on the opposite right and left sides of the 
dopants as imposed by the the bulk dimerisation |21|. If 
dopants are on the same sub-lattice, solitons are located 
on the same side of the dopants [22J and the effective 
exchange J eS (Ai) is ferromagnetic and decays rapidly to 
become negligible when Ai > 2£. The behavior of the 
GS local staggered magnetization is plotted on Fig. El for 
parameters corresponding to 4 typical behaviors. The 
key result here is the fact that the effective pairwise in- 
teraction is not frustrating (because of its sign alterna- 
tion with distance) although frustration is present in the 
microscopic underlying model. AF ordering is then ex- 
pected (at T = 0) as seen for a related system of coupled 
Spin-Peierls chains 0. 

We finish with a discussion on the relevance to ma- 
terials. Coexistence of AF ordering and SP order ob- 
served in Zinc-doped CuGe03 [lfj can be understood 
as an ordering of the induced 5=1/2 moments in the 
dimerised background @. We have shown that a 4-spin 
magnetic inter-chain interaction, similarly to the cou- 
pling to a 2D (or 3D) adiabatic lattice |£| but in contrast 
to a pairwise inter-chain coupling, leads to the formation 
of local magnetic moments under doping by inert atoms. 
Intermediate-range non-frustrated spin exchange interac- 
tions between these moments are expected to stabilize a 
low-tempcrature AF phase. Therefore, cyclic exchange 
is expected to play a key role in quasi-lD frustrated spin 
chains like L1V2O5 [2j. 

We gratefully acknowledge stimulating discussions 
with I. Affleck and E. Sorensen. 

Note added.-Anadditional small confinement cx J\j J\ 
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